Adiabatic compression of a trapped Fermi gas 
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We propose a method to reach conditions of high degener- 
acy in a trapped Fermi gas, based on the adiabatic transfer of 
atoms from a magnetic to a tighter optical trap. The transfor- 
mation yields a large increase of the Fermi energy, without a 
significant change of the temperature. The large enhancement 
of the central density emphasizes the role of the interactions 
and makes the system much closer to the BCS transition. An 
estimate of the time needed to achieve the conditions of adia- 
baticity is also given. 

The experimental realization of a highly degenerate 
atomic Fermi gas confined in traps is a task of primary 
importance, especially in view of the perspective of ap- 
proaching the BCS transition to the superfluid phase. 
The regime of quantum degeneracy has been already 
reached in a sample of potassium atoms where first 
signatures of Fermi statistics, like the deviation of the 
velocity distribution from a Boltzmann profile and the 
increase of the kinetic energy with respect to the clas- 
sical value, have been observed. The main difficulties 
in further lowering the temperature are due to the fact 
that the efficiency of the evaporative cooling process is 
strongly quenched In fact Fermi statistics inhibits 

collisional processes at low temperature both directly, by 
reducing the phase space available for collisions, and in- 
directly by lowering the density of the sample because of 
Pauli repulsion. Procedures to optimize the evaporative 
process have recently permitted to reach lower tempera- 
tures, of the order of 0.2 ~ 0.3 Tp where Tp is the Fermi 
temperature j|. 

The purpose of this work is to propose a method to 
reach conditions of high degeneracy, based on an adi- 
abatic compression of the gas. A similar method has 
proven quite successful in producing Bosc-Einstein con- 
densation in a reversible way, starting from a trapped 
Bose gas above the critical temperature Q]. The main 
point is that, by changing the shape of the confinement 
from a harmonic to a non harmonic trap, one can increase 
the degree of quantum degeneracy by keeping the entropy 
of the total system constant ||. At the same time, the 
process of thermalization is not drastically quenched and 
in typical experimental conditions takes place over times 
much shorter than the lifetime of the cloud. In the follow- 
ing we consider a gas occupying two different spin states, 
initially confined by a harmonic trap. We then switch on 
adiabatically a second tighter trap (see Fig. 1). Exper- 
imentally this can be realized using a magnetic trap for 



the first confinement and an optical trap for the second 
one. As a consequence of the adiabatic process a fraction 
of atoms will move from the magnetic to the optical trap. 
This can provide several important advantages: 

i) The gas in the optical trap becomes much more de- 
generate than the original one. In particular, if the num- 
ber of atoms transferred to the optical trap is a small 
fraction, the temperature will not change significantly 
with respect to the initial value, but the Fermi energy 
will increase, in a way proportional to the depth of the 
optical trap. 

ii) The gas in the optical trap is much denser due to 
the tighter confinement. This produces an increase of 
interaction effects and hence of the value of the critical 
temperature. Both effects i) and ii) favour the reachabil- 
ity of the BCS transition. 

iii) The Fermi energies of the two spin components, 
which initially were different because of the different mag- 
netic trapping, become closer in the optical trap, thereby 
favouring the mechanism of Cooper pairing. 

Another important advantage of the proposed method 
is that the velocity distributions of the atomic clouds oc- 
cupying the magnetic and optical traps can be measured 
separately. By releasing first the magnetic trap, one can 
measure the temperature of the sample. The effects of 
quantum degeneracy can then be investigated by mea- 
suring the velocity distribution of the gas confined in the 
optical trap. 

In the first part of the work we assume that thermo- 
dynamic equilibrium is ensured during each step of the 
adiabatic compression and we explore the properties of 
the new gas produced in the tight confinement by impos- 
ing entropy conservation. In the second part of the work 
we provide estimates for the times required to ensure 
adiabaticity and discuss possible scenarios where faster 
adiabatic trasformations take place out of thermal equi- 
librium. 

Let us consider a gas initially confined in a harmonic 
trap (hereafter called magnetic trap). We assume that 
the trapping frequencies and the number of atoms are 
the same for the two spin species. The Fermi energy is 
given by e F = huj mag (6N) 1 / 3 , where N is the number of 
atoms of each species and w mag is the geometrical average 
of the frequencies characterizing the magnetic trapping 
potential V ma g. We will consider systems lying initially 
in configurations of moderate degeneracy corresponding 
to ksTi — 0.2 — 0.5 e%, where Ti is the initial temperature 
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of the gas. For simplicity, we will assume that also the 
optical trap can be approximated by a harmonic potential 
V op t having a tighter frequency oj op t 3> w mag , and depth 
V opt (r = 0) = — U. In our model the trapping potential 
is then defined as V ext (r) = V opt (r) inside the optical 
trap (V op t < Vmag), and V ext (r) = V mag (r) outside (see 
Fig. 1). The maximum number of atoms that can be 
transferred in the optical trap is given by the value 
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As we will see, if we start from moderately low tempera- 
tures Eq. (|l|) provides an accurate estimate of the number 
of atoms which are actually transferred by the adiabatic 
process, provided N opt -C N. The relative number of 
atoms transferred in the optical trap is then given by the 
useful expression 
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Typical values that will be considered are Wmag/^opt — 
0.1 and U/e% = 5, corresponding to N opt /N ~ 10%. 
Since the number of transferred atoms is small one ex- 
pects that the final temperature Tf of the gas will not 
change significantly with respect to the initial value Ti. 
The final degree of degeneracy of the gas will be how- 
ever significantly higher, since the final Fermi energy is 
approximately given by ep ~ U + ep. At the same time 
the central density of the gas, that at small tempera- 
ture is equal to n(0) = (2me F /h 2 ) 3 ^ 2 /(67r 2 ) in the initial 
stage, will increase by the factor (ep /e°p) 3 ^ 2 in the op- 
tical trap. The increase of the Fermi energy and of the 
central density has an important effect on the value of 
the BCS temperature which is expected to behave, for 
negative scattering lengths, as 
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where pp — 7i(67r 2 ) 1 / 3 n(0) 1 / 3 is the Fermi momentum 
calculated in the center of the trap. The value of the 
BCS temperature can increase significantly with respect 
to its value in the magnetic trap. For example, by taking 
a r~j — 2000ao, a ratio U/e F = 5 and an initial central 
density n(0) ~ 10 12 cm -3 , we obtain an increase of the 
density by a factor ~15 and the BCS temperature (||) 
increases by the huge factor ~50, becoming comparable 
to the initial value of the Fermi temperature ep/ks- The 
above discussion suggests that the proposed adiabatic 
mechanism might provide conditions of high degeneracy, 
not far from the transition to the BCS phase. With such 
a denser gas also the effects of the mean field on the 
density profile may be significant. An estimate is given 
by the ratio [M Ei nt /Eh — 0.3pFa/h between the inter- 
action energy and the oscillator energy of a spherically 



symmetric trap. By using the values employed above one 
finds corrections of the order of 10 %. 

The high degeneracy realized in the gas is expected to 
show up in the velocity distribution and in the release 
energy. After completing the adiabatic transfer one can 
release the magnetic trap. Measuring the velocity distri- 
bution of these atoms then provides information on the 
temperature of the system which is expected to be close 
to the initial value. The atoms of the optical trap can 
be imaged in a second step. For them one predicts that 
the ratio Ekin/^BTf ~ 3U/8kBTf between the kinetic 
energy per particle and the thermal energy, should be 
enhanced in a significant way if U ksTf revealing the 
effects of quantum degeneracy. 

To confirm the scenario emerging from the above dis- 
cussion we have carried out a numerical calculation of the 
thermodynamic functions before and after the adiabatic 
transformation. The calculation is obtained by imposing 
that the initial and final configurations have the same 
entropy. This has been calculated using the semiclassical 
expression 
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where z = exp(/i//csT) is the gas fugacity and 



e(p,r) = ^- + V ext (r) 
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are the semiclassical particle energies. The results are 
presented in Figs. 2-5. In Fig. 2 we show the relative 
number of atoms in the optical trap as a function of the 
depth of the optical trap U/e F for two initial tempera- 
tures. For the chosen configuration with tiJmag/^opt = 0.1 
and final depth U = 5e° F , the optical trap can host about 
10 % of atoms confirming the analytic prediction (Q). 

In Fig. 3 the final temperature of the gas is plotted as 
a function oiU/e F . As already anticipated the final tem- 
perature of the gas does not change significantly from the 
original value, except for values of U/e F of the order of 
w pt/u)mag- Notice however that, due to the tight confine- 
ment, the temperature of the gas can become comparable 
to the optical oscillator temperature hio op t/kB if is not 
large enough. 

In Fig. 4 we show the kinetic energy per atom of 
the gas in units of the final temperature. The initial 
temperature is ksTi = 0.25e^. In the same figure, we 
also show the kinetic energy per particle of the gas oc- 
cupying separately the magnetic and the optical trap. 
These values are obtained by averaging separately the ki- 
netic energy over the particles occupying the optical trap 
(e(p,r) < 0), and the particles occupying the magnetic 
trap (e(p, r) > 0). The kinetic energy provides an impor- 
tant indicator of the quantum degeneracy of the gas. The 
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effect is very spectacular for the atoms of the optical trap, 
where at U = 5e F one finds Eki n /k B Tf ~ 7. By choosing 
a higher initial temperature the effect is less pronounced, 
but still large (E km /k B T f ~ 3 for k B Ti = 0.5e^ and 
U = 5e F ). Because of the high degeneracy the velocity 
distribution of the atoms in the optical trap deviates sig- 
nificantly from a Boltzmann distribution (Fig. 5). We 
have also calculated the central density as a function of 
the depth U. For U = 5e° F and k B T t = 0.25e° F , wc find 
an increase by a factor 15 with respect to the initial value. 

In the second part of the work we discuss the condi- 
tions needed to achieve the proposed adiabatic transfor- 
mation. Let us recall that there are several time scales in 
the problem. A first scale is fixed by the periods of the 
harmonic wells. These times (of order of 10~ 3 - 10 _1 sec) 
are expected to be shorter than the relaxation times due 
to collisions. The condition of reversibility requires that 
the relaxation time be much shorter than the time over 
which the optical trap is switched on. We have simulated 
the process of thermalization during the gradual increase 
of the depth of the optical trap by solving the quantum 
Boltzmann equation. We assume equal distribution func- 
tions for the two spin components and that the phase- 
space distribution of particles is a function only of the 
single-particle energies e(p, r) (ergodic assumption) 0,0. 
This yields the following equation to solve 

9/(ei) ma f 
P(ei) dt = ^3 / de z dt A p(e min ) 

x [/( e3 )/(e4)(l-/(e 1 ))(l-/(e 2 )) 

- /(e 1 )/(e 2 )(l-/(e 3 ))(l-/(e 4 ))] , (6) 

where p(e) is the density of states in the potential V ex t, 
which we model by p(e) — (e + U) 2 /2(huj pt) 3 if e < 
0, and p(e) = e 2 /2(huj ma g) 3 if e > 0. Furthermore, 
£min =min{ei, e2, £3, £4} is the minimum value of the 
four single-particle energies involved in the collisions and 
(2 = (3 + £4 — £i according to energy conservation. The 
cross-section for collisions between the two distinguish- 
able spin states is a = Ana 2 and is fixed by the s-wave 
scattering length a. In the numerical simulation the 
depth U of the optical trap is ramped up by keeping 
the ratio u) mag /uj op t constant, and at each step we let 
the gas thermalize in the new configuration. From Eq. 

, it turns out that the time scale of thermalization 
is fixed by = w mQff (|a|/a mag ) 2 iV 2 / 3 , where a mag = 
(h/mumag) 1 / 2 is the magnetic oscillator length. If wc 
choose uj mag /uj op t — 0.1, U = 5e F and k B Ti = 0.5e F , we 
find that the time required for the reversible transforma- 
tion is t rev ~ lOOr. For a configuration with TV = 10 6 , 
^mag — 2n x 100 Hz and \a\/a mag = 5x 10~ 3 , this corre- 
sponds to t rev ~ 1 sec. If one starts from the lower initial 
temperature k B Ti ~ 0.25e^ and keeps the other param- 
eters unchanged, the time t rev turns out to be about a 
factor two shorter. The origin of this behavior is due to 
the mechanism of atom exchange between the magnetic 



and the optical trap. Notice, however, that t rev depends 
rather strongly on the value of N and becomes longer if 
N decreases. 

The true shape of the optical trap differs from the har- 
monic potential employed above. Differences appear be- 
cause a dipole trap can be safely approximated by a har- 
monic potential only at its center and, in general, can 
host more atoms than the corresponding harmonic trap 
with the same depth U and frequency ratio Lo mag jui op t. 
This results in more heating and, as a consequence, in 
a higher final temperature Tf. For k B Ti — 0.25e F , we 
find that the fraction of atoms in the optical trap at 
U = 5e F is about 50% larger than the value obtained 
using the harmonic approximation, while the final tem- 
perature is only about 10% higher. The final degree of 
degeneracy of the gas in the optical trap is consequently 
slightly reduced, but the qualitative features discussed 
above remain unchanged. 

If the time of the transformation is longer than the in- 
verse oscillator frequency, but faster than the relaxation 
time (a relatively easy condition to realize experimen- 
tally), then the transformation will be adiabatic, in the 
sense of entropy conservation, but the system will not 
be in thermal equilibrium. Such a transformation cor- 
responds to an adiabatic transfer of the lowest N opt sin- 
gle particle states from the magnetic to the optical trap. 
The transformation keeps the corresponding occupation 
numbers unchanged. In this case, the atoms in the op- 
tical trap form a cold, out of equilibrium gas. With the 
initial condition k B Ti — 0.5 this out-of-equilibrium 
transformation transfers approximately the same frac- 
tion of atoms into the optical trap as with the fully re- 
versible transformation discussed above and the average 
kinetic energy of these atoms takes a similar value, re- 
vealing that the system is highly degenerate. By using 
the quantum Boltzmann equation (^J) we have found that 
the time needed for this configuration to relax to equi- 
librium is significantly shorter, the final temperature of 
the gas being only sligthly higher than the one obtained 
in the reversible transformation. 

In conclusion we have investigated the consequences of 
an adiabatic transfer of a gas of Fermions from a mag- 
netic to a tighter optical trap. We have seen that it 
is possible to reach configurations involving a significant 
fraction of atoms which, as a result of the transformation, 
will occupy a Fermi sea in conditions of high degeneracy. 
The large enhancement of the release energy should be 
easily observable through time of flight measurements. 
This method could make it possible to produce highly 
degenerate Fermi gases, which, in the case of gases inter- 
acting with negative scattering length, are close to the 
transition to the BCS phase. 
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FIG. 1. Schematic representation of the confining potential 
V ex t(r) in the presence of the magnetic and optical traps. 
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FIG. 2. Fraction of atoms in the optical trap as a function 
of U/e F for two initial temperatures: fcgTi = 0.56^ (dashed 
line), and ksTi — 0.25e% (dotted line). The solid line corre- 
sponds to the analytical estimate (0). 
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FIG. 3. Ratio between the final and initial temperature of 
the function of U /e% for the same initial temperatures 

as in Fig. 2. 
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FIG. 4. Kinetic energy per particle Ekin in units of the 
final temperature ksTf as a function of U/e° F . The initial 
temperature is ksTi — 0.25e^. The solid line corresponds 
to the kinetic energy per particle averaged over the entire 
system. Also shown are the kinetic energy per particle in the 
optical trap (dashed line), and in the magnetic trap (dotted 
line) . 
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FIG. 5. Momentum distribution of the atoms occupying 
the optical trap with U = §e° F (solid line). The initial value 
of the temperature is /cbTi = 0.25eJ-. The dashed line corre- 
sponds to a Boltzmann gas of the same number of atoms at 
the temperature T = T/. 
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